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November 9, 2024 at 06:48 


1. Intro. This is a modification of the program SSMCC, extending it by adding proof logging. We will 
use the VeriPB program as the external verifier. This program creates two extra output files: an OPB file 
containing a system of linear inequalities to model our MCC problem, and a PBP file which contains our 
proof logging. the input format is the same as the one defined in the program DLX3. The proof logging 
extensions were introduced by Filip Stappers in 2024. 

This program is another experiment in the use of so-called sparse-set data structures instead of the dancing 
links. It is written as if living on a planet where the sparse-set ideas are well known, but doubly linked links 
are almost unheard-of. 

I suggest that you read SSXCC, SSXCC-BINARY, DLX3 and SSMCC first. 

After this program finds all solutions, it normally prints their total number on stderr, together with 
statistics about how many nodes were in the search tree, and how many “updates” were made. The running 
time in “mems” is also reported, together with the approximate number of bytes needed for data storage. 
(An “update” is the removal of an option from the list of one its items.) One “mem” essentially means a 
memory access to a 64-bit word. The reported totals don’t include the time or space needed to parse the 
input or to format the output.) 


7##define o memst+ /* count one mem */ 

#define 00 mems += 2 /* count two mems */ 

#-define 000 mems += 3 /* count three mems */ 

##define subroutine.overhead mems += 4 

#define O "/" /* used for percent signs in format strings «/ 
#define mod % /* used for percent signs denoting remainder in C */ 
##define maz_stage 500 /* at most this many options in a solution */ 


#define maz_level 32000 /* at most this many levels in the search tree «/ 

#define maz-_cols 100000 /* at most this many items */ 

7##define maz_nodes 10000000 /* at most this many nonzero elements in the matrix */ 
7##define savesize 10000000 /* at most this many entries on savestack */ 

#define bufsize (9 * maz_cols + 3) /* a buffer big enough to hold all item names «/ 


#define show_basics 1 /* vubose code for basic stats; this is the default «/ 

#define show_choices 2 /* vbose code for backtrack logging */ 

#define show-_details 4 /* vbose code for further commentary */ 

##define show_record_weights 16 /* vbose code for first time a weight appears */ 
#define show_weight_bumps 32 /* vubose code to show new weights */ 

#define show_final_weights 64 /* vubose code to display weights at the end «/ 
#define show_profile 128 /* ubose code to show the search tree profile */ 

##define show_fullstate 256 /* vbose code for complete state reports */ 

#define show_tots 512 /* vbose code for reporting item totals at start */ 

#define show_warnings 1024 /* vubose code for reporting options without primaries */ 
##define show_maz_deg 2048 /* vbose code for reporting maximum branching degree */ 
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2. Here is the overall structure: 


#include <stdio.h> 

#include <stdlib.h> 

#include <string.h> 

#include <ctype.h> 

#include "gb_flip.h" 
typedef unsigned int uint; /* a convenient abbreviation */ 
typedef unsigned long long ullng; /* ditto «/ 


(Type definitions 8 ); 
(Global variables 3); 
(Subroutines 6); 
int main(int argc,char xargv[]){ register int c, cc, 1, j, k, p, pp, q, 7, 8, t, cur-choice, cur_node, 
best_itm, istage, score, bests, best_l, line , optnr; 
(Process the command line 4); 
(Input the item names 15); 
(Input the options 19); 
(Output our problem to the OPB file for proof logging 30); 
(Output header of the PBP file for proof logging 35); 
if (vbose & show_basics) (Report the successful completion of the input phase 27); 
if (vbose & show_tots) (Report the item totals 28); 
imems = mems,mems = 0; 
if (baditem) (Report an uncoverable item 25 ) 
else { 
if (randomizing) {Randomize the item list 29); 
(Solve the problem 44); 
} 
done: (Write contradiction line for proof logging 40); 
if (vbose & show_profile) (Print the profile 64); 
if (vbose & show_final.weights) { 
fprintf (stderr, "Final,weights:\n"); 
print_weights (); 


if (vbose & show_maz_deg) fprintf(stderr, "The maximum best_itm size,was,"O"d.\n", mazdeg); 
if (vbose & show_basics) { 
forintf (stderr, "Altogether,."O"1llu,solution"O"s,"O"llut+"O"1lu,mems,", count, 
count =1?"":"s",imems,mems); 
bytes = (itemlength + setlength) * sizeof (int) + last_node * sizeof 
(node) + 2 * maal « sizeof (int) + mazsaveptr * sizeof (threeints ); 
fprintf (stderr, "\"O"1lluvupdates, "O"1lubytes, "O"llu,nodes,", updates , bytes, nodes); 
fprintf (stderr, "\yccosty"O"11d%%.\n", mems ? (200 * cmems + mems)/(2 * mems) : 0); 
} 
if (sanity_checking) fprintf(stderr, "sanity_checking,was,on!\n"); 
(Close the files 5); 
} 
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3. You can control the amount of output, as well as certain properties of the algorithm, by specifying 


options on the command line: 


e ‘v(integer )’ enables or disables various kinds of verbose output on stderr, given by binary codes such as 


show-_choices ; 


e ‘m( integer )’ causes every mth solution to be output (the default is mO, which merely counts them); 


? 


e ‘s( integer ) 


causes the algorithm to randomize the initial list of items (thus providing some variety, 


although the solutions are by no means uniformly random); 

e ‘d( integer )’ sets delta, which causes periodic state reports on stderr after the algorithm has performed 
approximately delta mems since the previous report (default 10000000000); 

e ‘c( positive integer )’ limits the levels on which choices are shown during verbose tracing; 

e ‘C( positive integer )’ limits the levels on which choices are shown in the periodic state reports (default 10); 


‘1(nonnegative integer )’ gives a lower limit, relative to the maximum level so far achieved, to the levels 


on which choices are shown during verbose tracing; 
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Global variables 3) = 


‘t (positive integer )’ causes the program to stop after this many solutions have been found; 

‘T( integer )’ sets timeout (which causes abrupt termination if mems > timeout at the beginning of a level); 
‘w( float )’ is the initial increment dw added to an item’s weight (default 1.0); 

‘W( float )’ is the factor by which dw changes dynamically (default 1.0); 

‘S$ (filename )’ to output a “shape file” that encodes the search tree. 
‘P( filename )’ (required) name of the 


6 


“opb file” and “.pbp file” for proof logging; 


int random_seed = 0; /* seed for the random words of gb_rand */ 

int randomizing; /* has ‘s’ been specified? */ 

int vbose = show_basics + show_warnings; /* level of verbosity */ 

int spacing; /* solution k is output if k is a multiple of spacing */ 

int show_choices.maz = 1000000; /* above this level, show_choices is ignored */ 

int show_choices.gap = 1000000; /* below level mazl — show-choices.gap, show_details is ignored */ 


int show_levels.maz = 10; 


/* above this level, state reports stop */ 


int mac; /* maximum level actually reached */ 

int mazs; /* maximum stage actually reached */ 

int mazsaveptr ; /* maximum size of savestack */ 

char buf [bufsize]; /* input buffer «/ 

ullng count; /* solutions found so far */ 

ullng options; /* options seen so far */ 

ullng imems, mems, tmems, cmems; /* mem counts */ 
ullng updates; /* update counts */ 

ullng bytes; /* memory used by main data structures */ 

ullng nodes; /* total number of branch nodes initiated */ 


ullng thresh = 10000000000; 
ullng delta = 10000000000; 


/* report when mems exceeds this, if delta £0 */ 
/* report every delta or so mems */ 


ullng maxcount = *ffffffffffff{fffFf; /* stop after finding this many solutions */ 
ullng timeout = *1f£ffffffffffffff; /* give up after this many mems */ 
float w0 = 1.0, dw = 1.0, dwfactor = 1.0; /* initial weight, increment, and growth */ 


float mazwt = 1.0; /* largest weight seen so far */ 
FILE xshape_file; /* file for optional output of search tree shape */ 
char xshape_name; /* its name x*/ 


FILE xopb_file; 


/* file for OPB model of the problem we are solving */ 


char opb_name [255]; /* its name x/ 


FILE xpbp-_file; 


/* file for our proof log */ 


char pbp_name [255]; /* its name «/ 


int mazdeg; 


/* the largest branching degree seen so far */ 


int plLconstraints; /*x Current number of proof logging constraints written */ 


See also sections 9, 43, and 45. 
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This code is used in section 2. 


4. If an option appears more than once on the command line, the first appearance takes precedence. 
(Process the command line 4) = 
for (j = argc —1,k =0; 7; j——-) 
switch (argv [j][0]) { 


case ’v’: k |= (sscanf (arguv|j] + 1,""O"d", &ubose) — 1); break; 

case ’m’: k |= (sscanf (argv|[j] +1,""O"d", & spacing) — 1); break; 

case ’s’: k |= (sscanf (argu|j] + 1,""O"d", &random_seed) — 1), randomizing = 1; break; 
case ’d’: k |= (sscanf (arguv[j] +1,""O"11d", &delta) — 1), thresh = delta; break; 
case ’c’: k |= (sscanf (argu[j] + 1,""O"d", &show_choices.maxz ) — 1); break; 
case ’C’: k |= (sscanf (argu[j] + 1,""O"d", &show_levels.maz ) — 1); break; 

case ?1’: k |= (sscanf (argu[j] + 1,""O"d", &show_choices_gap ) — 1); break; 
case ’t’: k |= (sscanf (argu[j] + 1,""O"11d", &maxcount) — 1); break; 

case ’T’: k |= (sscanf (argv[j] + 1,""O"11d", &timeout) — 1); break; 

case ’w’: k |= (sscanf (arguv[j] +1,""O"£", &dw) — 1); break; 

case ’W’: k |= (sscanf (arguv[j] + 1,""O"£", &dwfactor) — 1); break; 

case ’S’: shape._name = argv|j] + 1, shape_file = fopen (shape_name, "w"); 


if (ashape_file ) 
forintf (stderr, "Sorry, ,Ijcan’ tyopen file,‘"O"s’ for writing! \n", shape._name); 
break; 
case ’P’: strcpy(opb_name, argv[j] + 1), strcat(opb_name,".opb"), opb_file = fopen(opb_name, "w"); 
if (sopb_file) fprintf (stderr, "Sorry, ,Ican’ tyopen file, ‘"O"s’ for writing! \n", opbname); 
strepy (pbp_name, argv |j] + 1), strcat (pbp_name,".pbp"), pbp_file = fopen(pbp_name, "w"); 
if (=pbp_file) fprintf (stderr, "Sorry, ,Ican’tyopen file, ‘"O"s’ for writing! \n", pbp_name); 
break; 
default: k = 1; /* unrecognized command-line option */ 
} 
if (k V sopb_file V apbp_file) { 
fprintf (stderr, "Usage: ,"O"s ,P<bar>,[v<n>] |, [m<n>] ,, [s<n>] ,[d<n>] "",[c<n>] ,[C<n>] , [1<n\ 
>] W[t<n>] | [T<n>] | [w<f>] |, [W<f>] |, [S<bar>] \<_f00.d1x\n", argv [0]); 
exit (—1); 
} 
if (randomizing) gb-init_rand (random-_seed ); 


This code is used in section 2. 


5. (Close the files 5) = 

if (shape_file) fclose (shape_file ); 
if (opb_file) fclose (opb_file); 

if (pbp_file) fclose (pbp_file); 


This code is used in section 2. 


6. Here’s a subroutine that I hope is never invoked (except maybe when I’m debugging). 


(Subroutines 6) = 
void confusion(char *m) 


{ 


forintf (stderr, ""O"s!\n",m); 


} 


See also sections 11, 12, 138, 14, 31, 41, 49, 55, 61, 62, and 63. 


This code is used in section 2. 


87 SSMCC-PROOF DATA STRUCTURES i) 


7. Data structures. Sparse-set data structures were introduced by Preston Briggs and Linda Torczon 
[ACM Letters on Programming Languages and Systems 2 (1993), 59-69], who realized that exercise 2.12 
in Aho, Hopcroft, and Ullman’s classic text The Design and Analysis of Computer Algorithms (Addison— 
Wesley, 1974) was much more than just a slick trick to avoid initializing an array. (Indeed, TAOCP exercise 
2.2.6—24 calls it the “sparse array trick.” ) 

The basic idea is amazingly simple, when specialized to the situations that we need to deal with: We can 
represent a subset S of the universe U = {2,%1,...,%n—1} by maintaining two n-element arrays p and gq, 
each of which is a permutation of {0,1,...,n — 1}, together with an integer s in the range 0 <s <n. In 
fact, p is the inverse of g; and s is the number of elements of $. The current value of the set S' is then simply 
{Lpo.+-+,Lp,_,}. (Notice that every s-element subset can be represented in s!(n — s)! ways.) 

It’s easy to test if x, € S, because that’s true if and only if gq, < s. It’s easy to insert a new element x, 
into S: Swap indices so that p, = k, q, = s, then increase s by 1. It’s easy to delete an element x; that 
belongs to S: Decrease s by 1, then swap indices so that p, = k and q, = s. And so on. 

Briggs and Torczon were interested in applications where s begins at zero and tends to remain small. 
In such cases, p and q need not be permutations: The values of ps, ps41, ---,; Pn—1 Can be garbage, and 
the values of gq, need be defined only when x, € S. (Such situations correspond to the treatment by Aho, 
Hopcroft, and Ullman, who started with an array full of garbage and used a sparse-set structure to remember 
the set of nongarbage cells.) Our applications are different: Each set begins equal to its intended universe, 
and gradually shrinks. In such cases, we might as well maintain inverse permutations. The basic operations 
go faster when we know in advance that we aren’t inserting an element that’s already present (nor deleting 
an element that isn’t). 

Many variations are possible. For example, p could be a permutation of {xo,21,...,2%n—1} instead of a 
permutation of {0,1,...,n —1}. The arrays that play the role of g in the following routines don’t have 
indices that are consecutive; they live inside of other structures. 
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8. This program has an array called item, with one entry for each item. The value of item[k] is an 
index x into a much larger array called set. The set of all options that involve the kth item appears in that 
array beginning at set[x]; and it continues for s consecutive entries, where s = size(a) is an abbreviation 
for set|a — 1]. If ttem[k] = x, we maintain the relation pos(x) = k, where pos(x) is an abbreviation for 
set[z — 2]. Thus item plays the role of array p, in a sparse-set data structure for the set of all currently 
active items; and pos plays the role of q. 

A primary item x also has a wt field, set[a — 5], initially 1. The weight is increased by dw whenever 
we backtrack because x cannot be covered. (Weights aren’t actually used in the present program; that will 
come in extensions to be written later. But it will be convenient to have space ready for them in our data 
structures, so that those extensions will be easy to write.) 

And finally, we have the bound and slack fields, set |x —6] and set|[a —7]. These are used in the same way 
as in DLX3 to keep track of the item’s multiplicities. 

Suppose the kth item x currently appears in s options. Those options are indices into nd, which is an 
array of “nodes.” Each node has three fields: itm, loc, and clr. Ifa <q < a+, let y = set|g]. This is 
essentially a pointer to a node, and we have nd[y].ttm = x, nd[y].loc = q. In other words, the sequential 
list of s elements that begins at « = item|k] in the set array is the sparse-set representation of the currently 
active options that contain the kth item. The clr field nd[y].clr contains x’s color for this option. The itm 
and clr fields remain constant, once we’ve initialized everything, but the loc fields will change. 

The given options are stored sequentially in the nd array, with one node per item, separated by “spacer” 
nodes. If y is the spacer node following an option with ¢ items, we have nd[y].itm = —t. If y is the spacer 
node preceding an option with t items, we have nd[y].loc = t. 

This probably sounds confusing, until you can see some code. Meanwhile, let’s take note of the invariant 
relations that hold whenever k, g, x, and y have appropriate values: 


pos(item|k]) =k; nd[set[q]].loc =q; item[pos(x)] =x; set|nd[y].loc] = y. 


(These are the analogs of the invariant relations p[q[k]] = q[p[k]] = & in the simple sparse-set scheme that 
we started with.) 

The set array contains also the item names. 

We count one mem for a simultaneous access to the i#m and loc fields of a node. We don’t count mems 
for accesses to the linenr field, since these are only relevant for proof logging. 


#define size(x) set[(x) —1].2 /* number of active options of the kth item, x */ 

#define pos(x) set|(x) — 2].i /* where that item is found in the item array, k */ 

#define Iname(x) set|(x) — 4] /* the first four bytes of 7’s name */ 

#define rname(x) set|(x) —3 /* the last four bytes of x’s name */ 

#define slack(x) set{(x) — 5].i /* if multiplicity [u..v], slack is equal to v-u and does not change «/ 
#define bound(x) set|(x) — 6]. /* residual capacity of this item */ 

#define wt(x) set[(x) — 7].f /* the current floating-point “weight” of x */ 


2 
a 


#tdefine primeztra 7 /* this many extra entries of set for each primary item */ 
7##define secondeztra 4 /* and this many for each secondary item */ 
7#define mazextra 7 /* maximum of primeztra and secondextra */ 


#tdefine ipropcount 6 /* the number of bytes used for each item in the input phase */ 
(Type definitions 8) = 
typedef struct node_struct { 


int itm; /* the item x corresponding to this node «/ 

int loc; /* where this node resides in x’s active set */ 

int clr; /* color associated with item x in this option, if any */ 

int linenr; /* line number from the input where this node comes from */ 
} node; 


typedef union { 
int 7; /* an integer (32 bits) «/ 
float f; /* a floating point value (fits in 4 bytes) */ 
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} tetrabyte; 


See also section 10. 


This code is used in section 2. 


9. (Global variables 3) += 


node nd[maz_nodes]; /* the master list of nodes */ 
int last_node; /* the first node in nd that’s not yet used x/ 
int item|maz_cols]; /* the master list of items */ 
int second = maz_cols; /* boundary between primary and secondary items */ 
int last_itm; /* items seen so far during input, plus 1 */ 
tetrabyte set|maz_nodes + mazextra * maz-_cols]|; /* active options for active items «/ 
int itemlength; /* number of elements used in item */ 
int setlength; /* number of elements used in set */ 
int active; /* current number of active items */ 
int oactive; /* value of active before swapping out current-choice items */ 
int baditem; /* an item with no options, plus 1 */ 
int osecond; /* setting of second just after initial input */ 
int force |maz_cols]; /* stack of items known to have size = bound — slack x/ 
int forced; /* the number of items on that stack «/ 
10. We're going to store string data (an item’s name) in the midst of the integer array set. So we’ve got 


to do 


some type coercion using low-level C-ness. 


(Type definitions 8) += 
typedef struct { 
int 1, 7; 
} twoints; 
typedef struct { 
int 1, s, 0; 
} threeints; 
typedef union { 
unsigned char str|8]; /* eight one-byte characters */ 
twoints Ir; /* two four-byte integers */ 
} stringbuf; 
stringbuf namebuf ; 


11. 


(Subroutines 6) += 


void print_item_name (int k, FILE x«stream) 


{ 


namebuf .lr.l = Iname(k), namebuf .Ir.r = rname(k); 
forintf (stream, "\"O" .8s8", namebuf .str); 
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12. An option is identified by the names of the items it contains. For proof logging we use the line number 
where the option appears in the input as the identifier of the option. Here is a routine that prints an option, 
given a pointer to any of its nodes. It also prints the position of the option in its item list. 


(Subroutines 6) += 
void print_option (int p, FILE *stream,int showpos ) 
{ 
register int k, q, 2; 
x = nd[p].itm; 
if (p > last.node Va <0) { 
forintf (stderr, "Illegal,option,"O"d!\n", p); 
return; 
} 
for (q=p;; ) { 
print_item_name(a, stream); 
if (nd[q].clr) fprintf (stream, ":"O"c", nd[q].clr); 
qt; 
x = nd[q].itm; 
if (x1 <0) { 
qt+= 2; 
fprintf (stream, "\,(Line,"O"d)", nd[q].linenr); 
x = nd[q].itm; 
} 
if (¢ =p) break; 
} 
k = nd{q}.loc; 
if (showpos > 0) fprintf (stream, "C"O"d of "O"d) \n",k — x +1, size(x)); 
else if (showpos = 0) fprintf (stream,"\n"); 
} 
void prow (int p) 
4 
print_option(p, stderr, 1); 


} 
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13. When I’m debugging, I might want to look at one of the current item lists. 
(Subroutines 6) += 
void print_itm (int c) 
{ 
register int p; 
if (c < primextra V c > setlength V pos(c) < 0V pos(c) > itemlength V item[pos(c)] 4c) { 
forintf (stderr, "Illegaljitem,"O"d!\n",c); 
return; 


forintf (stderr, "Item,("O"d)",c); 
print.item_name (c, stderr); 
if (c < second) { 
if (slack(c) V bound(c) £1) fprintf(stderr, "\("O"d,"O"d)", bound(c) — slack (c), bound (c)); 
forintf (stderr, "WC"O" dof y"O"d) , ylength"O"d, wweight,"O".1£:\n", pos(c) + 1, active, 
size(c), wt(c)); 


else if (pos(c) > active) 
fprintf (stderr, "\,(secondary,"O"d, jpurified) , ,length,"O"d:\n", pos(c) + 1, size(c)); 
else fprintf(stderr, "\(secondary,"O"d) , ylength,"O"d:\n", pos(c) + 1, size(c)); 
for (p=c; p<c+size(c); p++) prow(set[p].i); 
} 
void print_items() 


{ 


register int 7; 


for (i =0; i < ttemlength; i++) print itm (item|i]); 


} 
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14. Speaking of debugging, here’s a routine to check if redundant parts of our data structure have gone 


awry. 
7##define sanity_checking 0 /* set this to 1 if you suspect a bug */ 


(Subroutines 6) += 
void sanity() 


register int k, x, 7, 1, r, q, 99; 
int ok = 1; 
for (k =0; k < itemlength; k++) { 
x = item|k); 
if (pos (i) £k) { 
fprintf (stderr, "Bad. posyfield of item"); 
print.item_name (a, stderr); 
fprintf (stderr, "WC"O"dy!=,"O"d, y"O"d) !\n", k, pos (x), x); 
ok = 0; 
} 


for (4 =0; i < last_node; i++) { 
1 = ndli].itm,r = nd{i].loc; 
if (I< 0) { 
if (nd[it+r+1].itm # —-r) { 
forintf (stderr, "Bad spacerin,nodes,"O"d,"O"d!\n", 7,4 +r +1); 
ok = 0; 
} 
qq = 90; 
} else { 
if (l>r) fprintf(stderr ,"itm>loc in node,"O"d!\n", 2); 
else { 
if (set[r].1 41) { 
fprintf (stderr, "Bad locyfield for option,"O"dofjitem",r —1+ 1); 
print_item_name (1, stderr); 
fprintf (stderr, "\yinunode,."O"d! , set [r] .i="O"d\n", i, set[r].i); 
ok = 0; 


if (pos(l) < active) { 

if (r <1+ size(l)) q=+1; else q=-—1; /* in or out? */ 

if (q* qq <0) { 
forintf (stderr, "Flipped statusatoption,,"O"d of item",r —1+ 1); 
print_item_name (I, stderr); 
forintf (stderr, "yin node,"O"d! ,4q,qq="O"d,"O"d\n", 2, g, 9q); 
ok = 0; 

} 


qq = | 
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15. Inputting the matrix. Brute force is the rule in this part of the code, whose goal is to parse and 
store the input data and to check its validity. 
We use only ipropcount entries of set per item, while initially reading the item-name line. 
#define panic(m) 
{ fprintf(stderr, ""O"s!\n"O"d:4"O".99s\n",m, p, buf); exit (—666); } 
(Input the item names 15) = 
line = 0; while (1) { 
if (fgets (buf , bufsize, stdin)) break; 
line ++ ; 
if (0, buf [p = strlen(buf) — 1] 4 ?\n’) panic("Inputylineway,jtoo, long"); 
for (p = 0; 0, isspace (buf [p]); p++) ; 
if (buf |p] =’ 1? V abuf [p]) continue; /* bypass comment or blank line «/ 
last_itm = 1; 
break; } 
if (slast_itm) panic("No,items"); 
for (; 0, buf [p]; ) { 
o, namebuf .lr.l = namebuf .lr.r = 0; 
(Scan an item name, possibly prefixed by bounds 16); 
00, name (last_itm * ipropcount ) = namebuf .lr.l, rname (last_itm * ipropcount) = namebuf .Ir.r; 
0, slack (last_itm * tpropcount) = q — 1, bound (last_itm * ipropcount) = q; 
last_itm ++; 
if (last_itm > maz_cols) panic("Too,many,,items"); 
for (p += j +1; 0, isspace (buf [p]); p++) ; 
if (buf [pl =? 1?) { 
if (second 4 maz_cols) panic("Itemname,line,contains,,| twice"); 
second = last_itm; 
for (p++; o, isspace (buf [p]); p++) ; 
} 
} 


This code is used in section 2. 
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16. (Scan an item name, possibly prefixed by bounds 16) = 
if (second = maz_cols) istage =0; else istage = 2; 
start_name: for (j = 0; 7 < 8A (0, nisspace (buf [p + j])); JH) { 
if (buf [p+ j)=?2?) { 
if (istage) panic("Illegal,,‘:’yin,item mame"); 
(Convert the prefix to an integer, q 17); 
r=q,tstage =1; 
goto start_name; 
} else if (buf[p+ 7] =’1’) { 
if (istage > 1) panic("Illegal,‘ |’ ,in,item name"); 
(Convert the prefix to an integer, g 17); 
if (q=0) panic("Upper_bound,is,zero" ); 
if (istage =0) r=q; 
else if (r > q) panic("Lower bound jexceeds,upper,bound"); 
istage = 2; 
goto start_name; 


} 


0, namebuf .str{j] = buf [p + J]; 


switch (istage) { 

case 1: panic("Lower,bound_without_upper,bound"); 

case 0: g=r=1; 

case 2: break; 

} 

if (j =0) panic("Itemnamejempty"); 

if (7 = 8A 7isspace (buf [p+ j])) panic("Item name, ,too,long"); 
(Check for duplicate item name 18); 


This code is used in section 15. 


17. (Convert the prefix to an integer, g 17) = 
for (q=0, pp =p; pp <p+Jj; pp++) { 
if (buf [pp] < °0’ V buf [pp] > °9’) panic("Illegal,digit,,in, bound spec"); 
q=10*q+ buf [pp] — 70’; 


p= pp +1; 
while (j) namebuf .str[——j] = 0; 


This code is used in section 16. 


18. (Check for duplicate item name 18) = 
for (k = last_itm — 1; k; k--) { 
if (0, name(k * ipropcount) 4 namebuf .Ir.l) continue; 
if (rname(k * ipropcount) = namebuf .Ir.r) break; 
if (k) panic("Duplicate,item mame" ); 


This code is used in section 16. 
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19. (Input the options 19) = 
while (1) { 
if (fgets (buf , bufsize, stdin)) break; 
line + ; 
if (0, buf [p = strlen(buf) — 1] 4 ?\n’) panic("Optionline,,too,,long" ); 
for (p = 0; 0, isspace (buf [p]); p++) ; 
if (buf [p] =’ |? V abuf [p]) continue; /* bypass comment or blank line */ 
i = last_node; /*x remember the spacer at the left of this option */ 
for (pp = 0; buf [p]; ) { 
o, namebuf .Ir.l = namebuf .lr.r = 0; 
for (j =0; j < 8A (0, 7isspace (buf [p + j])) A buf [p+ j] 42:73 G++) 0, namebuf .str|j] = buf [p + J]; 
if (47) panic("Empty item name"); 
if (j =8A 7isspace (buf [p + j]) A buf [p+ 7] 4? +’) panic("Item mame, too,,long"); 
(Create a node for the item named in buf [p] 20); 
if (buf [p+ j] 4 ’?:’) 0, nd[last_node].clr = 0; 
else if (k > second) { 
if ((0, isspace (buf [p + 7 + 1])) V (0, misspace (buf [p + 7 + 2]))) 
panic("Color,must,be,a,single, character"); 
o, nd[last_node].clr = (unsigned char) buf [p+ 7 + 1); 
pt+=2; 
} else panic("Primaryitem must jbe,uncolored"); 
for (p += j +1; 0, isspace (buf [p]); p++) ; 


} 
if (spp) { 
if (vbose & show_warnings) fprintf (stderr, "Option ignored,,(no.primary,jitems) :,"O"s", buf); 
while (last_node >i) { 
(Remove last_node from its item list 21); 


last_node —-; 
} else { 
o, nd[i].loc = last_node — 1; /* complete the previous spacer */ 
last_node ++; /* create the next spacer */ 
if (last_node = maz_nodes) panic("Too,jmany nodes"); 
options ++; 


o, nd[last_node].itm = i+ 1 — last_node; 


Expand set 23); 
Adjust nd 24); 
(Make optionless items invisible 26 ); 


} 
\ (Initialize item 22); 
( 
( 


This code is used in section 2. 
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20. We temporarily use pos to recognize duplicate items in an option. 


(Create a node for the item named in buf [p] 20) = 
for (k = (last_itm — 1) * ipropcount; k > 0; k —= ipropcount) { 
if (0, name(k) 4 namebuf.ir.l) continue; 
if (rname(k) = namebuf .Ir.r) break; 
} 
if (4k) panic("Unknown,item, name"); 
if (0, pos(k) > 7%) panic("Duplicate,item,mamein,this,option"); 
last_node-++; 
if (last_node = maa_nodes) panic("Too many modes"); 
0, t = size(k); /* how many previous options have used this item? */ 
o, nd[last_node|.itm = k/ipropcount, nd |last_node|.loc = t, nd[last_node].linenr = line ; 
if ((k/ipropcount) < second) pp = 1; 
o, size(k) =t +1, pos(k) = last_node; 


This code is used in section 19. 


21. (Remove last_node from its item list 21) = 
0, k = nd[last_node].itm * ipropcount; 
00, size(k)——, pos(k) =i-—1; 


This code is used in section 19. 


22. (Initialize item 22) = 

active = itemlength = last_itm — 1; 

for (k = 0,7 = primextra; k < itemlength; k++) { 
oo, item|k] = j,j += (k +2 < second ? primeztra : secondextra) + size((k + 1) * tpropcount ); 
if (j < item|k] + ipropcount) j = item|k] + ipropcount; 

} 

setlength = j — ipropcount; /* a decent upper bound x/ 

if (second = maz-cols) osecond = active, second = j; 

else osecond = second — 1; 


This code is used in section 19. 


§20 


23. Going from high to low, we now move the item names and sizes to their final positions (leaving room 


for the pointers into nb). 


(Expand set 23) = 
for (; k; k—) { 

0,j = item|k — 1); 
if (k = second) second = J; /* second is now an index into set */ 
00, size(j) = size(k * tpropcount); 
0, pos(j) =k — 1; 
00, rname(j) = rname(k * ipropcount), Iname(j) = name (k * ipropcount ); 
oo, slack(j) = slack (k * ipropcount ), bound(j) = bound (k * ipropcount); 
if (k < osecond) { 


o, wt(j) = w0; 
if (size(j) < bound(j) — slack(j)) baditem = k; 
else if (size(j) =0) force [forced ++] = J; 

} else if (size(j) =0) force[forced ++] = J; 


i 


This code is used in section 19. 
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24. (Adjust nd 24) = 
for (k =1; k < last_node; k++) { 


if (0, nd[k].itm <0) continue; /* skip over a spacer */ 
0,j = item|nd|[k].itm — 1]; 
i= 7 + nd[k].loc; /* no mem charged because we just read nd|[k].itm */ 


o, nd|k].ttm = j, nd[k].loc = 4; 
0, set |i].i = k; 

} 
This code is used in section 19. 


25. (Report an uncoverable item 25) = 
{ 
if (vbose & show_choices) { 
k = item[baditem — 1); 
forintf (stderr, "Item"); 
print_item_name(k, stderr); 
fprintf (stderr, "jhasyfewer than,"O"d options! \n", bound(k) — slack (k)); 
} 
i 


This code is used in section 2. 


26. (Make optionless items invisible 26) = 
while (forced) { 
0,j = force|—— forced]; 
if (vbose & show_details) { 
forintf (stderr, "Deactivating,optionless,item"); 
print_item_name (j, stderr); 
forintf (stderr, "\n"); 


00,1 = item|[—- active], pp = pos(j); 
00, item|active] = j, item[pp] = 4; 
00, pos(j) = active, pos(i) = pp; 


} 


This code is used in section 19. 


27. The “number of entries” includes spacers (because DLX2 includes spacers in its reports). If you want 
to know the sum of the option lengths, just subtract the number of options. 
(Report the successful completion of the input phase 27) = 
forintf (stderr, "("O"11ldjoptions,,"O"d+"O"d items, .,"O"d entries,successfullyread) \n", 
options , osecond, itemlength — osecond, last_node); 


This code is used in section 2. 
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28. The item lengths after input are shown (on request). But there’s little use trying to show them after 
the process is done, since they are restored somewhat blindly. (Failures of the linked-list implementation in 
DLX2 could sometimes be detected by showing the final lengths; but that reasoning no longer applies.) 


(Report the item totals 28) = 
{ 
forintf (stderr, "Item totals:"); 
for (k =0; k < itemlength; k++) { 
if (k = second) fprintf (stderr,"|"); 
forintf (stderr, "\"O"d", size (item[k])); 
} 
fprintf (stderr, "\n"); 
} 


This code is used in section 2. 


29. (Randomize the item list 29) = 
for (k = active; k>1; ) { 
mems += 4,7 = gb_unif_rand(k); 
k--} 
00, 00,t = item|j], item[j] = item|[k], item|k] = t; 
00, pos(t) = k, pos (item[j]) = J; 
i 


This code is used in section 2. 
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30. Proof logging. We will convert the MCC problem to a pseudo-boolean problem, so it can be verified 
by VeriPB. In the pseudo-boolean problem, all variables are booleans (0 or 1). We create one variable for 
each option (named Ojineny), and one variable for each item in each option. Primary items are named 
Olinenr-variablename, secondary items also include their color. 

Then we add constraints for each option/item pair, so that items are selected if and only if the option is 
selected. We add the constraint that for each primary item, the sum of all its corresponding ’item-in-option’ 
variables is 1, to make sure every primary item is covered. Finally, we add a constraint for each pair of 
incompatible options due to secondary items. 

We keep count of the number of constraints that we output during proof logging. Linear equalities count 
double, since they are converted to two linear inequalities by VeriPB. 

No mems are counted during proof logging. 

For example, let’s consider again the example from DLX3: 


A simple example of color controls 
B2:31IC | XY 

B X:0 Y:0 

C X:1 Y¥:1 
X 
X 
Y 


QAQAWarr e— 


ee © 


The first option will be transformed in these 4 lines, to ensure the option and its variables are both either 
true or false: 


1 o38_A 1 ~03 = 1; 
1 0o3.B 1 “03 = 1; 
1 038_X_0 1 ~03 = 1; 
1 o38_Y_.0 1 ~03 = 1; 


To ensure item A will be covered, this line is added: 
1 o8_A 1 04_A = 1; 


Option C has multiplicities. We convert these to two constraints: 


1 04.C 1 05C 1 o7_C >= 2; 
-1 04.C -1 05C -1 o7_C >= -3; 


Finally, to ensure item X isn’t given two different colors by the options of lines 3 and 4, we add this 
constraint: 
1 ~o3_X_0 1 ~04X_1 >= 1; 


(Output our problem to the OPB file for proof logging 30) = 
plconstraints = 0; 
(Output constraints for options and items 32); 
(Output constraints to ensure each primary item is covered 33); 
(Output constraints to ensure secondary items are assigned only one color 34); 
if (opb_file) fclose (opb_file); 


This code is used in section 2. 
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31. (Subroutines 6) += 
void print_pLitemname(FILE xstream,int optnr,int k,int clr) 


namebuf .lr.l = Iname(k), namebuf .lr.r = rname(k); 
forintf (stream, "0o"O"a_"O" .8s", optnr, namebuf .str); 
if (k > second) { 

if (clr) fprintf (stream,"_"O"c", clr); 

else fprintf (stream,"_[]"); 


} 


32. (Output constraints for options and items 32) = 
for (1 = 0; i < last_node; i++) { 
k = nd[i].itm; 
if (k >0) { 
forint (opb_file, "1u"); 
print_plitemname (opb-_file, nd[2].linenr, k, nd{i].clr); 


forint (opb_file, "41470" O"dy=u13;\n", nd{i].linenr , namebuf .str, nd[i].linenr); 
pl.constraints += 2; 


} 
} 


This code is used in section 30. 


33. 


{ 


(Output constraints to ensure each primary item is covered 33) = 


int ww, ss, opt; 
for (k =0; k < active; k++) { 
ii = item|k]; 
if (ii > second) continue; 
for (s = ti, ss = s+ size(w#); s < ss; st+) { 
opt = set|s].2; 
printf (opb_file,"14."); 
print_pLitemname ( opb_file, nd[opt].linenr, ii, 0); 
fprintf (opb_file, "\."); 


if (slack (ii) =0) { 
fprintf (opb_file, "="O"d;\n", bound (ii)); 


else { 
fprintf (opb_file, ">="O"d;\n", bound (ii) — slack (ii)); 
for (s = ii, ss = +4 size(t); 6 < ss; st+) { 
opt = set[s].1; 
forintf (opb_file, "-1."); 
print_plitemname (opb-_file, nd[opt].linenr, ii, 0); 
forintf (opb-file, ".."); 


fprintf (opb_file, ">=,-"O"d;\n", bound (ii)); 


plconstraints += 2; 
} 
} 


This code is used in section 30. 
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34. For the secondary items, we add a condition in the OPB file for each pair of incompatible options. 


(Output constraints to ensure secondary items are assigned only one color 34) = 


{ 
int i, sl, s2, ss, optl, opt2; 
for (k =0; k < active; k++) { 
ii = item|k]; 
if (ii < second) continue; 
for (si = ti,ss = s1 + size(ti); 51 < ss; s1++) { 
optl = set[s1].i; 
for (s2 =s1+1; s2 < ss; s2++) { 
opt2 = set|s2].1; 
if (nd[opt1].clr 4 nd[opt2].clr V nd[opt1].clr =0V nd[opt2].clr =0) { 
fprintf (opb_file, "147"); 
print_pLitemname ( opb_file, nd[opt1 ].linenr, item|k], nd [opt |.clr); 
forintf (opb_file, "41u~"); 
print_pLitemname ( opb_file, nd[opt2].linenr, item[k], nd|[opt2|.clr ); 
forintf (opb_file, "4>=1;\n"); 
plconstraints + ; 


This code is used in section 30. 


35. We also need to write the necessary steps to the OPB file, so that VeriPB can verify our proof. This 
uses the routines below. 
(Output header of the PBP file for proof logging 35) = 

forint (pbp_file, "pseudo-Boolean, proof,,version,1.0\n"); 

forint (pbp_file, "£"O"d\n", pl_constraints ); 


This code is used in section 2. 


36. (Write a solution for proof logging 36) = 
sprintf (pbp_file,"vu"); 
int opt; 
for (k =0; k < stage; k++) { 
for (opt = choice [levelstage[k]]; nd[opt].itm > 0; opt——) ; 
forintf (pbp_file, "o"O"d,.", nd[+ opt].linenr); 
for (opt; nd[opt].itm > 0; opt++) { 
print_plitemname (pbp-file, nd[opt].linenr, nd[opt].itm, nd[opt].clr); 
Sprintf (pbp-file, "u"); 


} 
printf (pbp_file,"\n"); 
pl_constraints ++; 


This code is used in section 57. 
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37. When backtracking, we need to add a RUP constraint. 
(Write RUP constraint on failed include_option for proof logging 37) = 


{ 
fprintf (pbp_file, "rup.s"); 
fprintf (pbp_file, "147 0"O"d,.", nd[cur_choice].linenr); 
for (k = 0; k < stage; k++) { 
fprintf (pbp_file, "1,7 0"O"d", nd[choice [levelstage [k]]].linenr); 


fprintf (pbp_file, ">=,1;\n"); 
plconstraints ++; 


} 


This code is used in section 44. 


38. (Write RUP constraint on backtrack for proof logging 38) = 
{ 
if (stage >0) { 
forintf (pbp_file, "rupu"); 
for (k =0; k < stage; k++) { 
fprintf (pbp_file, "147 0"O"d\", nd[choice [levelstage [k]]].linenr ); 


forintf (pbp_file, ">=,1;\n"); 
pl_constraints ++; 


} 
} 


This code is used in section 44. 


39. When purging options that have no support, we also need to log a RUP constraint. 
(Write RUP constraint for proof logging when detected unsupported option 39) = 


i 
int optt, k; 
for (optt = opt; nd[optt].itm > 0; optt——) ; 
fprintf (pbp_file, "rupi1y~o"O"dy", nd [(+4 optt)|.linenr); 
for (k =0; k < rup_stage; k++) { 
for (optt = choice [levelstage [k]]; nd[optt].itm > 0; optt——) ; 
forintf (pbp_file, "1,7 0"O"da,", nd[(++ optt)].linenr); 


fprintf (pbp-file, ">=1u;\n"); 
pl_constraints ++; 


} 


40. When we finish our program, we need to notify our proof logging that we should have derived a 
contradiction. 


( Write contradiction line for proof logging 40) = 
printf (pbp_file, "cu-1\n"); 


This code is used in section 2. 


841 SSMCC-PROOF PROOF LOGGING 21 


41. For better performance of the verifier, it is useful to remove derived constraints from deeper levels of 
our backtrack tree, once they become reduntant for the further verification. For each level of our backtrack 
tree, we log when we enter this level. When we backtrack, we wipe all constraints derived on the lower levels, 
since these have become reduntant. 


(Subroutines 6) += 
void write_proof_log_level (int level) 


fprintf (pbp_file, "#4"O"d\n", level); 
void wipe_proof_log_level (int level) 


forintf (pbp_file, "w4"O"d\n", level); 
} 


22 BINARY BRANCHING VERSUS D-WAY BRANCHING SSMCC-PROOF §42 


42. Binary branching versus d-way branching. Nodes of the search tree in the previous program 
SSXCC, on which this one is based, are characterized by the name of a primary item 7 that hasn’t yet been 
covered. If that item currently appears in d options {01,...,0a}, node ¢ has d children, one for each choice 
of the option that will cover 7. 

The present program, however, makes 2-way branches, and its nodes are labeled with both an item 7 and 
an option o. The left child of node (i,0) represents the subproblem in which 7 is covered by 0, as before. 
But the right child represents the subproblem for which option o is removed but item 2 is still uncovered 
(unless d = 1, in which case there’s no right child). Thus our search tree is now rather like the binary tree 
that represents a general tree. (See The Art of Computer Programming, Section 2.3.2.) 

There usually is no good reason to do binary branching when we choose 7 so as to minimize d. On the 
right branch, 7 will have d — 1 remaining options; and no item 7’ will have fewer than d — 1. 

But this program is intended to provide the basis for other programs, which extend the branching heuristic 
by taking dynamic characteristics of the solution process into account. While exploring the left branch in 
such extensions, we might discover that a certain item 7’ is difficult to cover; hence we might prefer to branch 
on an option o’ that covers 7’, after rejecting o for item 2. 


43. We shall say that we’re in stage s when we’ve taken s left branches. We’ll also say, as usual, that we’re 
at level 1 when we’ve taken / branches altogether. 

Suppose, for instance, that we’re at level 5, having rejected 0; for 71, accepted og for i2, accepted o3 for is, 
rejected o4 for i4, and rejected os for is. Then we will have stage = 2, and choice{k] = o, for 0 < k < 5; 
here each o; is a node whose itm field is 7;. Also 


stagelevel 


stagelevel 
levelstage [0] = 
stagelevel 


I 


1, 
levelstage [1] = 2, 
levelstage [2] = 5. 
stagelevel 


0 

1 

2 
stagelevel [3] = 2 

4 

5 


stagelevel 


The option choice[k] has been accepted if and only if levelstage|stagelevel [k]] = k. 


(Global variables 3) += 
int stage; /* number of choices in current partial solution «/ 
int level; /* current depth in the search tree (which is binary) */ 
int choice |maz_level]; /* the option and item chosen on each level */ 
int deg|maz_level]; /* the number of options the item had at that time */ 
int levelstage|maz_stage]; /* the most recent level at each stage */ 
int stagelevel |maz_level]; /* the stage that corresponds to each level «/ 
ullng profile|maz_stage]; /* number of search tree nodes on each stage */ 
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44. The dancing. Our strategy for generating all exact covers will be to repeatedly choose an active 
primary item and to branch on the ways to reduce the possibilities for covering that item. And we explore 
all possibilities via depth-first search. 

The neat part of this algorithm is the way the sets are maintained. Depth-first search means last-in-first- 
out maintenance of data structures; and the sparse-set representations make it particularly easy to undo 
what we’ve done at deeper levels. 

The basic operation is “including an option.” That means (i) removing from the current subproblem all 
of the other options with which it conflicts, and (ii) considering all of its primary items to have their bounds 
decreased by 1. If this would make the bound of an item 0, we can make that item inactive. 


(Solve the problem 44) = 
{ 
level = stage = 0; 
forward: nodes++; 
write_proof_log_level (level); 
if (vbose & show_profile) profile|stage|++; 
if (sanity_checking) sanity(); 
(Maybe do a forced move 54); 
(Do special things if enough mems have accumulated 46 ); 
(Set best_itm to the best item for branching, and let score be its branching degree 56); 
if (forced) { 
0, best_itm = force |—— forced]; 
(Do a forced move and goto advance 60); 
} 
if (score = inf.size) (Visit a solution and goto backup 57); 
(Save the currently active items and their sizes and bounds 58); 
advance: oo, choice [level] = cur_choice = set|best_itm].i; 
o, deg [level] = score; 
if (minclude_option(cur-choice)) { 
(Write RUP constraint on failed include_option for proof logging 37); 
goto tryagain; 
} 
(Increase stage 47); (Increase level 48 ); 
goto forward; 
tryagain: if (score =1) goto prebackup; 
if (vbose & show_choices) fprintf (stderr, "Backtracking,,instage,"O"d_,(1)\n", stage); 
goto purgeit; 
prebackup: o, saveptr = saved [stage]; 
backup: write_proof_log_level (levelstage |stage — 1]); 
(Write RUP constraint on backtrack for proof logging 38 ); 
if (deg [level] A 1 A levelstage [stage — 1] #0) { 
wipe_proof_log_level (level ); 
} 
if (——stage <0) goto done; 
if (vbose & show_choices) fprintf (stderr, "Backtracking,tostage,"O"d_,(2) \n", stage); 
o, level = levelstage | stage]; 
purgeit: if (0, deg[level] =1) goto prebackup; 
(Restore the currently active items and their sizes and bounds 59); 
0, cur_choice = choice [level]; 
(Remove the option cur_choice 52); 
(Increase level 48 ); 
goto forward; 


} 
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This code is used in section 2. 


45. We save the sizes of active items on savestack, whose entries have two fields / and r, for an item and 
its size. This stack makes it easy to undo all deletions, by simply restoring the former sizes and bounds. 


(Global variables 3) += 
int level; /* number of choices in current partial solution */ 
int choice |maz_level]; /* the node chosen on each level */ 
int saved|maz_level + 1]; /* size of savestack on each level */ 
threeints savestack|savesize]; 
int saveptr; /* current size of savestack */ 
int tough_itm; /* item whose set of options has just become empty */ 


46. (Do special things if enough mems have accumulated 46) = 
if (delta \ (mems > thresh)) { 
thresh += delta; 
if (vbose & show_full_state) print_state( ); 
else print_progress(); 
} 
if (mems > timeout) { 
fprintf (stderr, "TIMEOUT! \n"); goto done; 


This code is used in section 44. 


47. (Increase stage 47) = 
if (++stage > mars) { 
if (stage > maz_stage) { 
fprintf (stderr, "Too many,stages ! \n"); 
exit (—40); 
} 


mazs = stage; 


} 


This code is used in section 44. 


48. (Increase level 48) = 
if (++level > maal) { 
if (level > maz_level) { 
forintf (stderr, "Too,manylevels!\n"); 
exit (—4); 


maal = level; 


oo, stagelevel [level] = stage, levelstage|stage] = level; 


This code is used in section 44. 


§49 


49. 
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The include_option routine extends the current partial solution, by hiding option opt. In addition, 
it will cover any primary items in opt if their bound after hiding opt becomes 0. The routine returns 0, 
however, if that would make some other primary item uncoverable. (In the latter case, tough_itm is set to 


the item that was problematic.) 


(Subroutines 6) += 
int include_option(int opt) 


{ 


} 


register int c, optp, nn, nnp, ss, i, it, p, pp, S; 
subroutine_overhead ; 
if (vbose & show_choices) { 
forintf (stderr, "S"O"d:", stage); 
print_option (opt, stderr, 1); 
} 
for (; 0, nd[opt — 1].ttm > 0; opt—) ; /* move to the beginning of the option */ 
for (; 0, (ti = nd[opt].itm) > 0; opt++) { 
pp = ndlopt].loc; /* where opt appears in ii’s set */ 
0, p = pos (it); /* where ii appears in item */ 
if (p > active) { 
if (4% > second) continue; /* secondary item has been purified */ 
confusion ("active"); /* primary item of an active option must be active */ 


(Cover or commit item ii, decrease bound of item ii if primary, potentially deactivate it, or return 


0 50); 
} 


return 1; 
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50. We need to remove the options that conflict with opt from the sets of their items. 


(Cover or commit item ii, decrease bound of item 7 if primary, potentially deactivate it, or return 0 50) = 
d 
if (ii < second) oo, bound(ii)—-; 
if (ii > second V bound(ii) =0) { 
0, 88 = size(it); 
if (4i < second) c=0; else 0,c = nd[opt].clr; 
for (s = ii +.ss—1; s> ti; s—) 
if (s # pp) { 
0, optp = set|s].2; 
if (c=0V (0, nd[optp].clr #c)) (Remove optp from its other sets, or return 0 51); 


0, p = pos (tt); /* note that pos(ii) might have changed x/ 
0, Wt = item|—— active]; 
00, item|active] = ii, ittem[p] = iit; 
00, pos (tt) = active, pos (tt) = p; 
} else { 
0, 88 = size(ti) — 1; 
if (00,858 < bound(i) — slack(i)) { 
if ((vbose & show_details) A level < show_choices_maz / level > maal — show_choices_gap) { 
forintf (stderr, "\jcan’ tcover"); 
print_item_name (item [ii], stderr); 
forintf (stderr ,"\n"); 


tough_itm = tt; 
forced = 0; 
return 0; /* abort the deletion, lest ii be wiped out */ 


if (ss =0) { /* Just deactivate item ii, no hiding needed x/ 
0, tt = item|—— active]; 
oo, item|[active] = i, item|p] = iii; 
00, pos (ti) = active, pos (iii) = p; 

} else { 
00, nnp = set|t# + ss].i, size(ti) = ss; 
oo, set [ti + ss].i = opt, set[pp].2 = nnp; 
oo, ndlopt].loc = ti + ss,nd[nnp|.loc = pp; 
updates ++; 

} 

} 
} 


This code is used in section 49. 
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51. At this point optp points to a node of an option that we want to remove from the current subproblem. 
We swap it out of the sets of all its items, except for the sets of inactive secondary items. (These have been 
purified, and we shouldn’t mess with their sets.) 


(Remove optp from its other sets, or return 0 51) = 


{ 
register int nn, i, iii, p, pp, ss, nnp; 
for (nn = optp; o,nd[nn — 1].itm > 0; nn——) ; /* move to beginning of the option */ 
for (; 0, (% = nd[nn].ttm) > 0; nn++) { 
p = nd[nn].loc; 
if (p > second A (0, pos(iit)) > active) continue; /* ti already purified */ 
0, $8 = size(t#i) — 1; 
if (p < second) { 
if (00, ss < bound(ii) — slack(wi)) { 
if ((vbose & show_details) A level < show-choices.max / level > maal — show-choices_gap) { 
forintf (stderr, "Can’ tycover"); 
print.item_name (ii, stderr); 
fprintf (stderr ,", ysize="O"d, jbound="O"d, yslack="O"d, ju="O"d\n", ss, bound (ii), 
slack (ii), bound (ii) — slack (it)); 
} 


tough_itm = wi; 
forced = 0; 
return 0; /* abort the deletion, lest ii be wiped out */ 


} 
if (ss =0) { 
0, Wt = item|—— active]; 
0, pp = pos (it); 
if (vbose & show_details) { 
fprintf (stderr, "Empty,option,list ,,deactivating"); 
print_item_name (it, stderr ); 
fprintf (stderr, "\n"); 
i 
00, item|active] = ti, item[pp| = iii; 
00, pos (tt) = active, pos (ti) = pp; 


} 


if (ss > 0) { 
0, nnp = set [ti + ss].i; 
0, size(wt) = 8s; 
oo, set [ti + ss]. = nn, set[p].c = nnp; 
oo, nd[nn].loc = it + ss,nd[nnp].loc = p; 
updates ++; 

i 

} 
} 


This code is used in section 50. 
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52. (Remove the option cur_choice 52) = 


{ 


register int 7, iii, ss, p, nnp; 


for ( ; 0, nd[cur_choice — 1].itm > 0; cur_choice——) ; /* move to beginning */ 
for ( ; 0, (4 = nd[cur_choice].itm) > 0; cur_choice++) { 
p = nd[cur_choice].loc; 
if (p > second A (0, pos(it)) > active) continue; /* tw inactive «/ 
0, 88 = size(ti) — 1; 
if (p < second) { 
if (00, 8s < bound(ii) — slack(wi)) { 
if ((vbose & show_details) A level < show choices.max / level > maal — show-choices_gap) { 
forintf (stderr, "\jcan’ tycover"): 
print_item_name (item |ii], stderr); 
forintf (stderr, "\n"); 
} 


goto prebackup; 


‘i 
if (ss =0) { 
0, Wt = item|—— active]; 
0, pp = pos (it); 
if (vbose & show_details) { 
fprintf (stderr, "Null_move, deactivating"); 
print_item_name (it, stderr ); 
fprintf (stderr, "\n"); 
} 
00, item|active] = ti, item[pp| = tii; 
00, pos (tt) = active, pos (ti) = pp; 
} 
i 
if (ss >0) { 
00, nnp = set[ii + ss].i, size(a) = ss; 
oo, set [ti + ss].i = cur_choice, set[p].i = nnp; 
oo, nd[cur_choice].loc = tt + ss,nd[nnp].loc = p; 
updates ++; 
} 
} 
} 


This code is used in section 44. 
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53. Ifa weight becomes dangerously large, we rescale all the weights. 

(That will happen only when dwfactor isn’t 1.0. Adding a constant eventually “converges”: For example, 
if the constant is 1, we have convergence to 2!" after 2!’ — 1 = 16777215 steps. If the constant dw is .250001, 
convergence to 8.38861e+06 occurs after 25165819 steps!) 

(Note: I threw in the parameters dw and dwfactor only to do experiments. My preliminary experiments 
didn’t turn up any noteworthy results. But I didn’t have time to do a careful study; hence there might be 
some settings that work unexpectedly well. The code for rescaling might be flaky, since it hasn’t been tested 
very thoroughly at all.) 

#define dangerous 1-10°7p 
#define wmin 1-107°°p 


(Increase the weight of toughitm 53) = 
cmems += 2, 00, wt(tough_itm) += dw; 
if (vbose & show_record_weights \ wt(toughitm) > maxwt) { 
maxwt = wt (tough_itm); 
fprintf (stderr, ""O"8.1£", mazrwt); 
print.item_name(tough_itm, stderr); 
fprintf (stderr, "\"O"11ld\n", nodes); 
} 
if (vbose & show_weight_bumps) { 
print_item_name(tough_itm, stderr); 
fprintf (stderr, "\Wwty"O".1£\n", wt (toughitm)); 
} 
dw *= dwfactor ; 
if (wt(tough_itm) > dangerous) { 
register int k; 
register float t; 
tmems = mems; 
for (k =0; k < itemlength; k++) 
if (0, item[k] < second) { 
o,t = wt(item[k]) «1 -10~7°p; 
0, wt (item|k]) = (t < wmin ? wmin : t); 


dw x= 1-107; 

if (dw < wmin) dw = wmin; 
w0 *=1-1072%p; 

if (w0 < wmin) w0 = wmin; 
cmems += mems — tmems; 


} 


54. At level 0, the force stack contains primary items that had no options. Their lower bound was 0, so 
they should simply not appear. 
(Maybe do a forced move 54) = 
while (forced) { 
0, best_itm = force|—— forced]; 
if (0, pos(best_itm) > active) continue; 
(Do a forced move and goto advance 60); 


r 


This code is used in section 44. 
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55. (Subroutines 6) += 
void print_weights (void) 
{ 
register int k; 
for (k =0; k < itemlength; k++) 
if (item|k] < second A wt(item[k]) 4 w0) { 
print_item_name (item |k], stderr ); 
fprintf (stderr ,"wty"O" .1£\n", wt (item|[k])); 
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56. The “best item” is considered to be an item that minimizes the branching degree. If there are several 
candidates, we choose the leftmost — unless we’re randomizing, in which case we select one of them at 
random. 

Consider an item that has four options {w, x,y,z}, and suppose its bound is 3. If the slack is zero, we’ve 
got to choose either w or xz, so the branching degree is 2. But if slack = 1, we have three choices, w or x or 
y; if slack = 2, there are four choices; and if slack > 3, there are five, including the “null” choice. 

In general, the branching degree turns out to be 1+ s — 6+ 1, where / is the length of the item, b is the 
current bound, and s is the minimum of b and the slack. This formula gives degree < 0 if and only if / is 
too small to satisfy the item constraint; in such cases we will backtrack immediately. (It would have been 
possible to detect this condition early, before updating all the data structures and increasing level. But that 
would make the downdating process much more difficult and error-prone. Therefore I wait to discover such 
anomalies until item-choosing time.) 

Let’s assign the score 1+ s—b+1 to each item. If two items have the same score, I prefer the one with 
smaller s, because slack items are less constrained. If two items with the same s have the same score, I 
(counterintuitively) prefer the one with larger b (hence larger 1), because that tends to reduce the size of the 
final search tree. 

Consider, for instance, the following example taken from MDANCE: If we want to choose 2 options 
from 4 in one item, and 3 options from 5 in another, where all slacks are zero, and if the items are 
otherwise independent, it turns out that the number of nodes per level if we choose the smaller item first is 
(1,3,6,6-3,6-6,6-10). But if we choose the larger item first it is (1,3,6,10,10-3,10-6), which is smaller 
in the middle levels. 

Notice that a secondary item is active if and only if it has not been purified (that is, if and only if it hasn’t 
yet appeared in a chosen option). 


#define inf_size *7f£ffffff 
(Set best_itm to the best item for branching, and let score be its branching degree 56) = 


{ 
score = inf_size, tmems = mems; 
if ((vbose & show_details) A level < show_choices_max / level > maal — show_choices_gap ) 
forintf (stderr, "Level"O"d:", level); 
for (k =0; k < active; k++) 
if (0, item[k] < second) { 
0, 8 = slack (item|[k]); 
if (0,s > bound (item|k])) s = bound (item |[k]); 
0, t = size(item[k]) + s — bound (item[k]) + 1; 
if ((vbose & show_details) A level < show_choices_max / level > maal — show_choices_gap) { 
print_item_name (item |k], stderr); 
if (bound (item[k]) A1Vs 40) { 
fprintf (stderr,"("O"d:"O"d,"O"d)", bound (item [k]) — s, bound (item |k]), t); 
} else fprintf(stderr ,"("O"d)",t); 


} 
if (¢=1) 
for (t = bound (item|k]) — slack (item[k]); 1 > 0; i-—) 0, force[forced ++] = item[k]; 
else if (t < score \(t < score V (s < best_s A\(s < best_s V (size (item|[k]) > best_l A (size (item|[k]) > 
best_l V (item |[k] < best_itm))))))) 
score = t, best_itm = item|k], best_s = s, best_l = size (item|k]); 


if ((vbose & show_details) A level < show_choices_maz / level > maal — show_choices_gap) { 
if (score = inf_size) fprintf(stderr, "jsolution\n"); 
else if (forced) { 
fprintf (stderr, "found,,"O"d forced: ", forced); 
for (4 =0; i < forced; i++) print_item_name (force |i], stderr); 
fprintf (stderr, "\n"); 
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} else { 
forintf (stderr, "joranching,,on"); 
print.item_name (best_itm, stderr ); 
fprintf (stderr, "("O"d)\n", score); 
i 
} 
if (score > maxdeg A score < inf_size \ aforced) maxdeg = score; 
if (shape_file) { 
if (score = inf_size) fprintf (shape_file, "sol\n"); 
else { 
fprintf (shape_file,"""O"d", score); 
print_item_name (best_itm, shape-_file ); 
fprintf (shape_file,"\n"); 


fflush (shape_file); 


cmems += mems — tmems; 


} 


This code is used in section 44. 


57. (Visit a solution and goto backup 57) = 
ft 

count ++; 

(Write a solution for proof logging 36); 

if (spacing A (count mod spacing =0)) { 
printf (""O"11d:\n", count); 
for (k =0; k < stage; k++) print_option (choice [levelstage [k]], stdout , 0); 
fflush (stdout); 


if (count > mazrcount) goto done; 
goto backup; 


} 


This code is used in section 44. 


58. (Save the currently active items and their sizes and bounds 58) = 
o, saved [stage] = saveptr; 
if (saveptr + active > maxsaveptr) { 
if (saveptr + active > savesize) { 
forintf (stderr, "Stack, overflow, (savesize="O"d) !\n", savesize); 
exit (—5); 
} 
mazsaveptr = saveptr + active; 
} 
for (p = 0; p< active; pt+) { 
mems += 4, savestack|saveptr + p].l = item |p], savestack|saveptr + p].s = size(item|[p]); 
if (item|[p] < second) o, savestack|saveptr + p|.b = bound (item [p]); 
} 
saveptr += active; 


This code is used in section 44. 
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59. (Restore the currently active items and their sizes and bounds 59) = 
o, active = saveptr — saved |stage]; 
saveptr = saved|stage]; 
for (p = 0; p< active; pt) { 
000, size (savestack |saveptr + p].l) = savestack |saveptr + p].s; 
if (savestack|saveptr + p].l < second) o, bound(savestack |saveptr + p].l) = savestack|saveptr + p].b; 


} 


This code is used in section 44. 


60. <A forced move occurs when best_itm has bound —slack remaining options. In this case we can streamline 
the computation, because there’s no need to save the current active sizes and bounds. (They won’t be looked 
at.) 


(Do a forced move and goto advance 60) = 
{ 
if ((vbose & show_choices) A level < show_choices_max) fprintf(stderr," (forcing) \n"); 
0, saved|stage] = saveptr; /* nothing is placed on savestack */ 
score = 1; 
goto advance; 


} 


This code is used in sections 44 and 54. 


61. (Subroutines 6) += 
void print_savestack (int start,int stop) 
{ 
int k; 
for (k = start; k < stop; k++) { 
print_item_name (savestack [k].1, stderr); 
forintf (stderr, "("O"d) ,u"O"ay"O"d\n", savestack [k].1, savestack[k].s, savestack [k].b); 
} 
} 


62. (Subroutines 6) += 
void print_state (void ) 
{ 
register int [, s; 
forintf (stderr, "Currentystate,(level,"O"d) :\n", level); 
for (l= 0; | < level; I++) { 
if (levelstage |stagelevel [l|] 41) fprintf(stderr ,"~"); 
print_option (choice [I], stderr, —1); 
fprintf (stderr, "\(ofy"O"d) \n", deg [I]); 
if (1 > show_levels.max) { 
forintf (stderr, "...\n"); 
break; 
} 
I 


forintf (stderr, "\y"O"11ldjsolutions,,"O"1ldmems, ,and,max,jlevel,"O"d_so jfar.\n", count, 
mems, mazl); 
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63. During a long run, it’s helpful to have some way to measure progress. The following routine prints a 
string that indicates roughly where we are in the search tree. The string consists of node degrees, preceded 
by ‘~’ if the node wasn’t the current node in its stage (that is, if the node represents an option that has 
already been fully explored — “we’ve been there done that”). 

Following that string, a fractional estimate of total progress is computed, based on the naive assumption 
that the search tree has a uniform branching structure. If the tree consists of a single node, this estimate is .5. 
Otherwise, if the first choice is the kth choice in stage 0 and has degree d, the estimate is (k —1)/(d+k-1) 
plus 1/(d+k-—1) times the recursively evaluated estimate for the kth subtree. (This estimate might obviously 
be very misleading, in some cases, but at least it tends to grow monotonically.) 


(Subroutines 6) += 
void print_progress (void) 
{ 
register int J, Il, k, d, c, p, ds =0; 
register double /, jd; 
fprintf (stderr, "\after,"O"1lldmems:."O"1lldsols,",mems, count); 
for (f = 0.0, fd = 1.0,1 =0; I < level; I++) { 
if (I < show_levels_maz ) 
fprintf (stderr ,"\"O"s"O"d", levelstage [stagelevel [l]] =1? "" :"~", deg[l}); 
if (levelstage|stagelevel [l]] =1) { 
for (k = 1,d= deg[l], 11 =1—1; ll > 0A stagelevel [ll] = stagelevel[l]; k++,d++,ll—) ; 
fd x=d, f += (k —1)/fd; /* choice | is treated like k of d «/ 


if (1 > show_levels.max \ ads) ds = 1, fprintf (stderr,"..."); 


fprintf (stderr, "\"O" .5£\n", f +0.5/fd); 
} 


64. (Print the profile 64) = 


fprintf (stderr, "Profile:\n"); 
for (k =0; k < mazs; k++) fprintf(stderr ,""O"3d:"O"1ld\n", k, profile [k]); 
} 


This code is used in section 2. 
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65. Index. 

active: 9, 13, 14, 22, 26, 29, 33, 34, 49, 50, 51, 
52, 54, 56, 58, 59. 

advance: 44, 60. 


argc: 2, 4. 
argu: 2, 4. 
b: 10. 


backup: 44, 57. 

baditem: 2, 9, 23, 25. 

bestitm: 2, 44, 54, 56, 60. 

best_l: 2, 56. 

best_s: 2, 56. 

bound: 8, 9, 138, 15, 23, 25, 33, 50, 51, 52, 
56, 58, 59. 

buf: 3, 15, 16, 17, 19. 

bufsize: 1, 3, 15, 19. 


bytes: 2, 3. 
ce: 2, 13, 49, 63. 
Oe 22: 


choice: 36, 37, 38, 39, 48, 44, 45, 57, 62. 
clr: 8, 12, 19, 31, 32, 34, 36, 50. 
cmems: 2, 3, 53, 56. 

confusion: 6, 49. 

count: 2, 3, 57, 62, 63. 

cur.choice: 2, 37, 44, 52. 

cur_node: 2. 

d: 63. 

dangerous: 53. 

deg: 43, 44, 62, 63. 


delta: 3, 4, 46. 
done: 2, 44, 46, 57. 
ds: 63. 


dw: 3, 4, 8, 53. 
dwfactor: 3, 4, 53. 
exit: 4, 15, 47, 48, 58. 


f: 8, 63. 
fclose: 5, 30. 
fd: 63. 
fflush: 56, 57. 
fgets: 15, 19. 
fopen: 4. 


force: 9, 23, 26, 44, 54, 56. 

forced: 9, 23, 26, 44, 50, 51, 54, 56. 

forward: 44. 

forintf: 2, 4, 6, 11, 12, 18, 14, 15, 19, 25, 26, 
27, 28, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 
41, 44, 46, 47, 48, 49, 50, 51, 52, 53, 55, 56, 
58, 60, 61, 62, 63, 64. 

gb_init_rand: 4. 

gb-rand: 3. 

gb_unif_rand: 29. 

a: 2, 8, 13, 14. 
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wt: 33, 34, 49, 50, 51, 52. 

wt: 49, 50, 51, 52. 

imems: 2, 3. 

include_option: 44, 49. 

inf_size: 44, 56. 

tpropcount: 8, 15, 18, 20, 21, 22, 23. 

isspace: 15, 16, 19. 

istage: 2, 16. 

item: 8, 9, 18, 14, 22, 23, 24, 25, 26, 28, 29, 33, 
34, 49, 50, 51, 52, 53, 55, 56, 58. 

itemlength: 2, 9, 13, 14, 22, 27, 28, 53, 55. 

itm: 8, 12, 14, 19, 20, 21, 24, 32, 36, 39, 43, 
49, 51, 52. 

ee 

I: 10, 14, 62, 63. 

last.itm: 9, 15, 18, 20, 22. 

last.node: 2,9, 12, 14, 19, 20, 21, 24, 27, 32. 

level: 41, 43, 44, 45, 48, 50, 51, 52, 56, 60, 62, 63. 

levelstage: 36, 37, 38, 39, 43, 44, 48, 57, 62, 63. 

linenr: 8, 12, 20, 32, 33, 34, 36, 37, 38, 39. 

il: 68. 

Iname: 8, 11, 15, 18, 20, 23, 31. 

loc: 8, 12, 14, 19, 20, 24, 49, 50, 51, 52. 

ir: 10, 11, 15, 18, 19, 20, 31. 

m: 6. 

main: 2. 

max.cols: 1, 9, 15, 16, 22. 

maz_level: 1, 43, 45, 48. 

max.nodes: 1, 9, 19, 20. 

maz_stage: 1, 43, 47. 

maxcount: 3, 4, 57. 

maxdeg: 2, 3, 56. 

mazextra: 8, 9. 

maal: 2, 3, 48, 50, 51, 52, 56, 62. 

mars: 3, 47, 64. 

mazsaveptr: 2, 3, 58. 

marwt: 3, 53. 

mems: 1, 2, 3, 29, 46, 53, 56, 58, 62, 63. 

mod: 1, 57. 

namebuf: 10, 11, 15, 16, 17, 18, 19, 20, 31, 32. 

nb: 23. 

nd: 8, 9, 12, 14, 19, 20, 21, 24, 32, 33, 34, 36, 
37, 38, 39, 49, 50, 51, 52. 


nn: 49, 51. 
nnp: 49, 50, 51, 52. 
node: 2, 8, 9. 


node_struct: 8. 
nodes: 2, 3, 44, 53. 
O: 1 

o: 1. 
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oactive: 9. 


ok: 14. 


oo: 1, 15, 21, 22, 23, 26, 29, 44, 48, 50, 51, 52, 53. 


ooo: 1, 59. 

opb_file: 3, 4, 5, 30, 32, 33, 34. 
opb.name: 3, 4. 

opt: 33, 36, 39, 49, 50. 
options: 3, 19, 27. 


optnr: 2, 31. 
optp: 49, 50, 51. 
optt: 39 

opti: 34 

opt2: 34 


osecond: 9, 22, 23, 27. 

p: 2, 12, 13, 49, 51, 52, 63. 

panic: 15, 16, 17, 18, 19, 20. 

pbp_file: 3, 4, 5, 35, 36, 37, 38, 39, 40, 41. 


pbp_name: 3, 4. 


pLconstraints: 3, 30, 32, 33, 34, 35, 36, 37, 38, 39. 
pos: 8, 13, 14, 20, 21, 23, 26, 29, 49, 50, 51, 52, 54. 


pp: 2,17, 19, 20, 26, 49, 50, 51, 52. 

prebackup: 44, 52. 

primestra: 8, 13, 22. 

printitem.name: 11, 12, 13, 14, 25, 26, 50, 51, 
52, 53, 55, 56, 61. 

printitems: 13. 

printitm: 13. 

print_option: 12, 49, 57, 62. 

print_pLitemname: 31, 32, 33, 34, 36. 

print_progress: 46, 63. 

print_savestack: 61. 

print_state: 46, 62. 

print.weights: 2, 55. 

printf: 57. 

profile: 43, 44, 64. 

prow: 12, 13. 

purgeit: AA. 

q: 2, 12, 14. 

qq: 14. 

r: 2, 10, 14. 

random seed: 4, 

randomizing: 3, 4. 

rname: 8, 11, 15, 18, 20, 23, 31. 

rup_stage: 39. 

s: 2, 10, 49, 62. 

sanity: 14, 44. 

sanity_checking: 2, 14, 44. 

saved: 44, 45, 58, 59, 60. 

saveptr: 44, 45, 58, 59, 60. 

savesize: 1, 45, 58. 

savestack: 1, 3, 45, 58, 59, 60, 61. 

score: 2, 44, 56, 60. 


3, 
2, 
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second: 9, 13, 15, 16, 19, 20, 22, 23, 28, 31, 33, 
34, 49, 50, 51, 52, 53, 55, 56, 58, 59. 

secondertra: 8, 22. 

set: 8, 9, 10, 18, 14, 15, 23, 24, 33, 34, 44, 
50, 51, 52. 


setlength: 2, 9, 13, 22. 
shape_file: 3, 4, 5, 56. 
shape.name: 3, 4. 
show_basics: 1, 2, 3. 


show_choices: 1, 3, 25, 44, 49, 60. 

show_choices.gap: 3, 4, 50, 51, 52, 56. 

show_choices_maz: 3, 4, 50, 51, 52, 56, 60. 

show_details: 1, 3, 26, 50, 51, 52, 56. 

show_final_weights: 1, 2. 

show-_full_state: 1, 46. 

show_levels.maz: 3, 4, 62, 63. 

show_maaz_deg: 1, 2. 

show-profile: 1, 2, 44. 

show_record_weights: 1, 53. 

show_tots: 1, 2. 

show_warnings: 1, 3, 19. 

show_weight_bumps: 1, 53. 

showpos: 12. 

size: 8, 9, 12, 18, 14, 20, 21, 22, 23, 28, 33, 34, 
50, 51, 52, 56, 58, 59. 

slack: 8, 9, 13, 15, 23, 25, 33, 50, 51, 52, 56. 

spacing: 3, 4, 57. 

ss: 33, 34, 49, 50, 51, 52. 

sscanf: A. 

stage: 36, 37, 38, 43, 44, 47, 48, 49, 57, 58, 59, 60. 

stagelevel: 43, 48, 62, 63. 

start: 61. 

startname: 16. 

stderr: 1, 2, 3, 4, 6, 12, 18, 14, 15, 19, 25, 26, 
27, 28, 44, 46, 47, 48, 49, 50, 51, 52, 53, 55, 
56, 58, 60, 61, 62, 63, 64. 


stdin: 15, 19. 

stdout: 57. 

stop: 61. 

str: 10, 11, 16, 17, 19, 31, 32. 
streat: A. 

strcpy: A. 


stream: 11, 12, 31. 
stringbuf: 10. 

strlen: 15, 19. 
subroutine_overhead: 1, 49. 
sl: 34. 

82: 34. 

t: 2, 53. 

tetrabyte: 8, 9. 
threeints: 2, 10, 45. 
thresh: 3, 4, 46. 
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timeout: 3, 4, 46. 

tmems: 3, 53, 56. 

toughitm: 45, 49, 50, 51, 53. 

tryagain: 44. 

twoints: 10. 

uint: 2. 

ullng: 2, 3, 48. 

updates: 2, 3, 50, 51, 52. 

vbose: 1, 2, 3, 4, 19, 25, 26, 44, 46, 49, 50, 
51, 52, 53, 56, 60. 

wipe_proof_log_level: Al, 44. 

wmin: 53. 

write_proof_log_level: 41, 44. 

wt: 8, 13, 23, 53, 55. 

w0: 3, 23, 53, 55. 

vw: 12, 14. 
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(Adjust nd 24) Used in section 19. 

(Check for duplicate item name 18) Used in section 16. 

(Close the files 5) Used in section 2. 

(Convert the prefix to an integer, q¢ 17) Used in section 16. 

(Cover or commit item ii, decrease bound of item %# if primary, potentially deactivate it, or return 0 50) 
Used in section 49. 

(Create a node for the item named in buf [p] 20) Used in section 19. 

(Do a forced move and goto advance 60) Used in sections 44 and 54. 

(Do special things if enough mems have accumulated 46) Used in section 44. 

(Expand set 23) Used in section 19. 

(Global variables 3, 9, 43, 45) Used in section 2. 

(Increase the weight of tough_itm 53) 

(Increase level 48 ) Used in section 44. 

(Increase stage 47) Used in section 44. 

(Initialize item 22) Used in section 19. 

(Input the item names 15) Used in section 2. 

(Input the options 19) Used in section 2. 

(Make optionless items invisible 26) Used in section 19. 

(Maybe do a forced move 54) Used in section 44. 

(Output constraints for options and items 32) Used in section 30. 

(Output constraints to ensure each primary item is covered 33) Used in section 30. 

(Output constraints to ensure secondary items are assigned only one color 34) Used in section 30. 

(Output header of the PBP file for proof logging 35) Used in section 2. 

(Output our problem to the OPB file for proof logging 30) Used in section 2. 

(Print the profile 64) Used in section 2. 

(Process the command line 4) Used in section 2. 

(Randomize the item list 29) Used in section 2. 

(Remove the option cur_choice 52) Used in section 44. 

(Remove last_node from its item list 21) Used in section 19. 

(Remove optp from its other sets, or return 0 51) Used in section 50. 

(Report an uncoverable item 25) Used in section 2. 

(Report the item totals 28) Used in section 2. 

(Report the successful completion of the input phase 27) Used in section 2. 

(Restore the currently active items and their sizes and bounds 59) Used in section 44. 

(Save the currently active items and their sizes and bounds 58) Used in section 44. 

(Scan an item name, possibly prefixed by bounds 16) Used in section 15. 

(Set best_itm to the best item for branching, and let score be its branching degree 56) Used in section 44. 

(Solve the problem 44) Used in section 2. 

(Subroutines 6, 11, 12, 13, 14, 31, 41, 49, 55, 61, 62, 63) Used in section 2. 

(Type definitions 8,10) Used in section 2. 

(Visit a solution and goto backup 57) Used in section 44. 

(Write RUP constraint for proof logging when detected unsupported option 39 ) 

(Write RUP constraint on backtrack for proof logging 38) Used in section 44. 

(Write RUP constraint on failed include_option for proof logging 37) Used in section 44. 

(Write a solution for proof logging 36) Used in section 57. 

(Write contradiction line for proof logging 40) Used in section 2. 
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